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PERMimVITY  PROBE  MODELLING 


1.  INTRODUCTION 


The  use  of  the  open-ended  probe  for  measurements  requires  a  semi-infinite  sample  to  be 
placed  in  contact  with  it  To  define,  in  practical  terms,  the  extent  of  the  sample  we  have  recently 
developed  a  finite  element  technique  to  Dxxlel  ti»  fringing  fields  of  open-ended  lines  and  their 
perturbation  by  a  material.  The  current  technique  accommodates  lossless  dielectric;  its 
extension  to  lossy  materials  will  be  briefly  outlined. 

Numerous  ptqiers  have  reported  the  use  of  open-ended  coaxial  lines  for  dielectric 
measurements  at  microwave  frequencies  (1-4).  The  technique  is  based  on  the  study  of  the 
fringing  admittance  of  the  open-ended  line  and  its  alteratitxi  by  a  lossy  dielectric.  The  system  is 
nondestructive;  its  advantages  are  of  a  practical  nature  mainly  associated  with  the  simple 
handling  of  samples.  The  principal  shoncoming  of  tiie  technique  is  diat  a  full  analysis  of  the 
end-of  line-admittaiKe  requires  extensive  numerical  calculation  and  does  not  lend  itself  to  real¬ 
time  measurements  (5)  while  simple  models  tiiat  assume  the  probe  to  be  nonradiating  and  only 
capacitive  in  nature  fail  progressively  at  high  frequency  and  fOT  lossy  materials.  In  this  report 
we  introduce  2  approximate  admittance  models  in  which  tiie  radiating  element  of  the  probe  is 
explicitly  treated  together  with  inductive  and  caqracitive  components.  The  models  will  be 
expressed  in  a  polynomial  form  suitable  for  a  numoical  solution  so  that  the  complex 
permittivity  of  the  sample  in  contact  with  the  probe  is  calculated  from  the  measured  values  of 
the  admittance.  A  comparative  study  of  the  performance  of  all  3  models  is  included.  Results 
show  that  the  2  new  nKxkls,  derived  from  the  same  set  of  assumptions,  are  not  equal  in 
perfcxmance.  The  better  of  the  2  models  is  used  in  conjunction  with  a  3-mm  probe  and  a 
Hewlett  Packard  (HP)  8720  network  analyzer  for  dielectric  measurements  on  standard  liquids. 
The  results  show  that  the  technique  is  suitable  for  measurements  on  lossy  materials 
to  20  GHz.  Its  accuracy  and  ptecision  are  assessed  by  comparing  the  measurements  with 
equivalent  values  from  the  literature. 
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2.  COAXIAL  PROBES  MODELLING 


A  second  order  polynomial  finite  element  technique  was  developed  in  1988  in  our  research 
group  (6)  to  model  the  fringing  capacitance  of  coaxial  probes.  The  implementation  of  the 
technique  requires  that  the  area  of  the  solution  be  divided  into  a  mesh  of  triangular  elements. 
The  potential  of  both  the  inner  and  outer  conductors  is  known  aitd  assumed  constant.  The 
potential  in  each  triangular  element  is  approximated  by  a  function  the  coefficient  of  which  is 
related  to  the  spatial  coordinates  of  the  element  Minimisation  of  the  potential  energy  of  the 
system  yields  an  approximation  of  the  potential  distribution  that  satisfies  Laplace's  equation. 

For  the  purpose  of  this  applicaticxi  a  coaxial  air  line  of  iiuier  and  outer  radii  a  and  b 
respectively  is  represented  as  in  Figure  1. 


Neumann  Boundary 


Figure  1.  Region  in  which  the  solution  is  implemented. 


The  technique  was  extensively  improved  in  the  course  of  the  current  project  In  particular, 
an  automatic  mesh  generator  was  developed  and  the  solution  used  to  plot  equipotential  lines  in 
the  area  under  study. 


2.1  Implementation  of  the  Solution 


The  region  of  the  solution  is  divided  into  triangular  elements.  The  mesh  is  partially 
bounded  by  Dirichlet  and  Neumann  boundaries  Figure  1. 

The  inner  and  outer  conductors  axe  held  at  1  and  0  volts  respectively. 

In  each  element  the  relative  permittivity  is  known  and  the  potential  <I>  is  approximated  by 
a  trial  function  with  constant  coefficients,  such  that 


(1) 


where  <l>j  is  the  potential  at  vertex  i  and  04  is  a  linear  expression  in  r  and  z,  the  vertex 
coordinates. 

In  terms  of  these  elements  the  potential  energy  stored  in  the  region  of  the  solution  is 
prc^rortiorud  to 


F  = 


(V4>)^  rdrdz 


(2) 


where  r^  is  the  radial  component  of  each  element  The  finite  element  method  was  implemented 
to  the  second  order  so  that  each  triangle  is  defined  by  6  nodes. 

Combining  Equations  1  and  2  we  obtain 


i  y  r 

F  =  27c  Z  EkZZi**  I  Vo,Vaj4»jrdrdz 

“I  H 


(3) 
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By  minimizing  Equation  3  with  respect  to  O  a  set  of  equations  is  foimed  which  can  readily 
be  solved  for  the  potential  at  each  node. 


Equipotential  lines  are  then  drawn  to  give  a  visual  representadon  of  the  fringing  field  of  the 
probe  and  its  perturbation  by  a  dielectric  medium. 


To  test  the  validity  of  the  rrxxlel  the  fringing  capacitance  of  the  probe  may  be  calculated 
using  the  solution  obtained  from  the  tixxiel 


C  = 


e«(V<|))2rdj.dz 

y2J., 


(4) 


(5) 


A  comparison  between  the  calculated  and  measured  value  of  the  capacitance  validates  this 
approach  to  the  problem. 


2.2  Automatic  Mesh  Generator 


At  the  onset  of  this  projea  the  creation  of  a  triangular  mesh  to  implement  the  finite  element 
technique  was  done  by  hand  before  each  calculation.  This  procedure  is  time  consuming  and 
prone  to  error  particularly  when  a  finer  mesh  is  required  to  improve  the  resolution  of  the 
model.  A  brief  or  otherwise  unsuitable  mesh  design  could  lead  to  inaccurate  results.  To 
overcome  these  difficulties  and  speed  up  the  process,  an  automatic  mesh  generator  was 
develqred. 

A  brief  mesh,  referred  to  as  the  basis  mesh,  is  made  in  the  region  of  solution  and  the 
potential  at  each  node  is  calculated  using  the  Fmite  Element  approximation.  The  basis  mesh 
needs  refining  particularly  in  the  regions  of  high  field  gradient  This  refining  is  accomplished 
by  comparing  the  average  potentials  of  adjacent  elements.  If  the  difference  is  greater  than  a 
predefined  limit,  a  node  in  both  of  the  2  elements  is  specified  using  the  angle  bisector  method. 
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Connecting  this  node  to  the  3  vertices,  the  conesponding  triangular  element  is  subdivided  into 
3  triangles  which  are  then  assigned  with  new  element  numbers.  After  all  the  elements  in  the 
problem  are  con:q)ared,  the  calculation  is  carried  out  again  in  the  new  mesh.  These  procedures 
are  continued  until  the  differences  in  potential  between  each  adjacent  elements  are  less  than  the 
set  limit  A  brief  mesh  with  10  elements  is  shown  in  Hgure  2  and  the  subdivision  of  elements 
3  and  8  is  shown  as  an  example. 


Figure  2.  The  subdivision  of  elements  3  and  8  by  the  angle  bisector  method. 


The  basis  mesh  can  also  be  segmented  by  another  technique,  namely  the  fractal 
segmentation.  In  this  method,  3  nodes  are  defined  at  the  midpoints  of  each  side  of  a  triangular 
element  Connecting  these  nodes  together,  4  equivalent  triangular  elements  are  formed  which 
are  similar  to  the  original  element  The  method  is  illustrated  in  Figure  3.  However,  unlike  the 
former  technique,  this  segmentation  must  be  applied  to  all  elements  in  the  system  at  the  same 
time  to  satisfy  the  boundary  conditions  of  adjacent  elements.  Hence  after  applying  the 
segmentation  once,  the  number  of  elements  in  the  new  mesh  is  four  times  that  of  the  original 
one. 
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2 


•  original  node 

*  additional  node 


Figure  3.  A  second  order  element  is  subdivided  into  4  equivalent  elements  by 
the  fractal  segmentation  technique. 


2.3  Equipotential  Contour  Plot 


Results  are  given  in  Figures  4*7.  Figure  4  demonstrates  the  potential  due  to  an  open- 
ended  air  line.  The  same  technique  was  used  to  generate  the  equipotential  lines  of  a  measuring 
probe  based  on  a  Teflon-filled  Vine  (Fig.  5),  an  air  line  tenninated  by  an  impedance  matched 
dielectric  window  (Fig.  6)  and  a  conical  probe  of  our  own  design  (Rg.  7).  This  latter  probe 
(UK  patent  application  882681.1)  was  designed  to  provide  deeper  penetration  and  to  increase 
the  sensitivity  of  measurements  at  frequencies  below  100  MHz. 
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Figure  4.  A  0.5  mm  inner  conduaor  radius,  air  line  radiating  into  a  dielectric  of  £=200. 
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Figure  6.  A  probe  with  internal  discontinuity  (i)  basic  configuration  and 

dimensions,  (ii)  the  mesh  used  in  region  (a),  (iii)  the  equipotendal 
plot  in  region  (b). 


9 


Figure  7.  A  conical  probe  (i)  basic  configuration  an  i  dimensions,  (ii)  the  mesh 
used  in  region  (a),  (iii)  the  equipotential  plot  in  region  (b). 


2.4  Fringing  Fields  in  Lossy  Media 


In  Figures  4-7  the  medium  in  contact  with  the  probe  was  assumed  real  (e'=constant,  e"  = 
0).  The  extension  of  the  technique  to  lossy  media  (e*  ^  0)  falls  outside  the  scope  of  this 

project  but  will  be  briefly  oudined.  If  the  electric  field  E  is  expressed  as  the  gradient  of  a  scaler 
potential 


E  =  - V<t> 


(6) 


Then  the  case  of  a  harmonic  field  Ampere’s  Law  will  be 


V  xH  =  -{<y+je*co)V(t> 


(7) 


where  H  is  the  magnetic  field,  a  the  conductivity ,  £  the  permittivity  and  to  the  radial  frequency. 
Taking  the  divergence  of  the  above  equation  gives 


(a  +  je*®)  V^((»  =  o 


(8) 


which  can  be  solved  using  our  existing  finite  element  technique  in  conjunction  with  the  newly 
developed  automatic  mesh  generator.  This  extension  to  our  programme  would  enable  the  study 
of  field  penetration  in  lossy  homogeneous  and  layered  media.  This  extension  would  also  be 
particularly  useful  in  defining  the  electric  field  penetration  while  performing  in  vivo 
measurements. 
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3.  DEVELOPMENT  OF  AN  ADMITTANCE  MODEL 


3.1  Theoretical  Background 


When  vector  network  analyzers  are  used  characterise  the  admittance  of  open-ended 
coaxial  probe^it  is  inhetendy  assumed  that  their  characteristics  at  a  given  frequency  and  under 
steady  state  conditions  are  indistinguishable  from  the  characteristics  of  a  linear  network  (La,  a 
combination  (series  or  shunt)  of  a  resistor  and  an  inductor  or  capacitor). 

In  the  case  of  an  open-ended  probe  the  appropriate  configuration  is  capacitive  susceptance 
B  in  parallel  with  a  conductance  G.  In  this  case,  it  is  move  convenient  to  express  its  equivalent 
impedance  in  terms  of  its  inverse,  the  admittance  Y,  thus 


Y  =  G+jB 


(9) 


where  G  and  B  ate  a  factor  of  the  propagation  constant  k  and  the  radial  dimensions  of  the  line. 
By  definition  the  complex  reflection  coefficient  of  the  probe  T  is  given  in  terms  of  this 
admittance  as 


r  =  (l-Y/Yo)/(l-FY/Yo)  ao) 

Yo  (=  1/Zo)  is  the  intrinsic  admittance  of  the  line.  The  complex  reflection  coefficient  F  is 
readily  measurable  using  the  experimental  setup.  The  emphasis  of  the  problem  at  this  stage  is 
the  relationship  of  the  admittance  Y  to  the  complex  permittivity  e*  of  the  sample  with  which  it 

is  in  contact.  This  contact  is  the  relationship  that  will  enable  e*  to  be  calculated  from  measured 
values  of  r  or  Y. 


The  relationship  between  F  and  e*  is  rather  complex  and  a  rigorous  analytical  expression  is 

not  available.  Several  authors  have  obtained  approximate  solutions  to  this  problem.  The  nxist 
complete  analysis  is  that  of  Mosig  et  al.(5),  who  calculated  the  admittance  of  the  open-ended 
line  by  matching  the  electromagnetic  fields  at  the  plane  of  the  aperture  incorporating  the 
fundamental  transverse  electromagnetic  (TEM)  and  higher  order  modes.  This  approach  while 
limited  by  certain  assumptions  is  the  most  accurate  model  available  for  an  end  of  line 
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admittance.  However,  this  analysis  cannot  be  perfomied  on  a  microcomputer  for 
measurements  in  real  time.  Three  alternative  models  have  been  studied. 


3.1.1  Equivalent  circuit  model  (Model  1) 


The  approach  adopted  in  the  first  instance  is  to  consider  Y  in  teims  of  its  lumped  equivalent 
circuit  parameters  The  susceptance  B  is  capacitive  made  up  of  internal  and  external  fringing 
fields  (Fig.  8): 


Y=G+jo)(Cf  +  Co(£*-jE”)) 


(11) 


where  (Zf  and  Cq  are  the  components  of  the  fringing  field,  to  is  the  angular  frequency  (rad/sec) 
and  e*  =  e’  -  je”  is  the  complex  permittivity  of  the  sample. 


Figure  8.  Model  I  approximation  of  the  coaxial  probe. 


13 


At  the  limit  of  low  frequencies  and  Cf  and  Cq  are  constants  that  depend  on  the 
dimensions  of  the  probe. 


Assuming  the  probes  are  used  at  frequencies  such  that  the  wavelength  of  the  propagating 
wave  is  significantly  larger  than  the  radial  dimensions  of  the  probe,  a  quasi-steady  state  would 
prevail.  From  these  equations  we  would  have: 

r  =  (1  -  jaJZo(Cf  +  Co  e*))  /  (1  +  £*))  (12) 


and 


e*  =  (1-  n  /  (jt«)CoZo(l+  O)  -  CfCo  (13) 

with  Cf  significantly  smaller  than  Cq.  Neglecting  the  contribution  of  Cf  to  the  permittivity 
would  introduce  an  errcn'  in  e*.  This  error  has  been  shown  to  be  small  compared  to  the  value 

of  e*  of  most  biological  samples  at  radio  and  microwave  frequencies.  This  led  some  authors  to 

simplifying  the  equivalent  circuit  of  the  open-ended  probe  to  a  single  shunt  capacitance  Cq  (1). 
In  this  case  the  reflection  coefficient  becomes 


r  =  (l-jwCoZo  £*)/(!+ jwCpZoE*)  (14) 


Implementation  of  Model  I 


Our  approach  to  the  problem  at  this  stage  is  to  consider  the  general  case  (Equation  11)  and 
assume  the  following: 

•  G  and  Cf  contribute  to  the  admittance  of  the  probe. 

•  Their  contribution  is  significantly  less  than  that  of  the  external  fringing  capacitance  C^, 

•  Their  contribution  may  exhibit  a  certain  degree  of  frequency  dependence  and  to  a  lesser 
extent  a  dependence  on  e*,  and 

•  Cq  may  also  exhibit  some  dependence  on  frequency  and  t* . 
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These  assumptions  lead  us  to  express  the  measured  leflecdon  coefficient  as: 


rn,=  r(aj.Co.e*)  +  F(co.Co.e*,Cf.G) 


(15) 


Under  certain  conditions,  notably  when  the  wavelength  of  the  source  is  much  larger  than 
the  radial  dimensions  of  the  probe,  the  correction  term  F  is  small  and  of  the  OTder  of  magnitude 
of  the  experimental  uncertainty.  It  is,  however,  good  experimental  practice  to  eliminate  it  from 
the  measured  parameter.  The  correaed  reflection  coefficient  F  may  then  be  solved  for  e*. 

Our  approach  to  this  correction  is  to  use  the  response  of  a  reference  sample  and  a  difference 
technique  (4)  to  normalise  the  unknown  response. 


3.1.2  Integral  expression  of  the  admittance  (Model  2) 

In  the  cylindrical  coordinate  system  (p,  0,  z),  Figure  9,  assuming  the  radial  electric  field 
intensity  over  the  aperture  Ep  (p,0)  is  inversely  proportional  to  p,  an  expression  fOT  aperture 
admittance  can  be  obtained  (7,8). 


YL=j  - 


/.Wo’'cos<>'^^^(«>'dP'dP 


(16) 


where : 


k^  =  (i>2poe* 

kj  =  a>2poeceo 

r2=p2+p2-2pp'cos<l'' 
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Yo  =  2rt 


(17) 


The  dielectric  window  is  assumed  lossless  of  pennittivity  e^,  e*  is  the  complex  permittivity 
of  a  semi-infinite  sample. 


t  2 


e* 


Figure  9.  Schematic  representation  of  the  coaxial  probe  in  terms  of  the  parameters 
of  Model  2  and  Model  3. 


Misra  (7)  has  shown  that  when  the  radial  dimensions  of  the  probe  are  much  smaller  than 
the  propagating  wavelength  Equation  (16)  can  be  approximated  by  the  first  terms  of  the  series 

expansion  for  the  exponential  term  and  used  to  calculate  e*.  numerically,  from  measured  values 
of  Y.  This  approximation  was  found  suitable  for  the  measurement  of  low  permittivity  materials. 
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Implementation  of  Model  2 
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'0 


C7= 


720 


1 

¥  ■ 


cos  0'[p^+p’^-2pp'cos  <t>’l^^d(t)'dp'dp 


^  .repo'^^,(b^-a6)(b^-a2) 

*  3360  '■  3 


+|<b'‘-a4)2] 


C9= 


40320 


cos  <|>’[p^+p’^-2pp'cos  ()>’]""d(|>’dp'dp 


The  second  fenn  of  the  series  expansion  goes  to  zero  on  integration  over  Since  the 
integrand  of  C,  is  singular  at  p  =  p*  for  <)>'  =  0,  it  is  calculated  analytically  and  can  be  expressed 
as 


C,  =  4(a+b)[E(m)  - 1] 


(19) 


where  E(m)  is  the  complete  elliptic  integral  of  the  second  kind,  which  can  be  approximated  by  a 
polynomial  (12),  and 


n,=^ah_^ 
(a  +  b)2 


(20) 


Since  Q,  where  i  =  0, ...,  10,  are  dependent  only  on  the  physical  dimensions  of  the 
aperture,  they  can  be  evaluated  fcH-  a  given  coaxial  line.  The  triple  integrals  are  performed 

numerically  by  Simpson's  three  points  rule.  The  intervals  [a,b]  and  [o,n]  are  subdivided  into 
2N  (-10)  subintervals  and  Simpson's  rule  is  applied  on  each  pair  of  consecutive  subintervals. 
The  maximum  tolerance  in  the  calculation  of  E(ro)  was  set  to  be  10’^®.  The  number  of  terms  in 
the  series  in  Equation  (18)  was  varied  between  3  and  20  to  allow  a  check  for  convergence  and 
accuracy  in  the  model.  The  values  of  Cj,  where  i  =  0, ...,  20,  were  calculated  for  a  50  £2  3^nm 
OD  probe.  These  values  were  then  substituted  into  equation  (18)  and  the  probe  admittance  Y 
was  evaluated  from  1 30  MHz  to  20  GHz.  As  the  number  of  Cj  terms  considered  was  increased 
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the  calculated  value  of  the  admittance  converged  to  a  constant  value  when  nine  terms  in  the 
series  expansion  were  included. 

Implicit  in  this  model  is  that  only  the  lowest  transmission  mode,  the  TEM  mode,  is 
propagated.  Although  this  mode  is  the  dominant  mode,  higher  order  modes  will  be  generated  at 
the  probe  sample  interface. 

In  a  measurement  situation  a  Newton-Raphson  technique  is  used  to  calculate  the  complex 
permittivity  of  the  sample  fiom  the  expeiimentaUy  measured  values  of  Y  at  each  frequency  of 
interest 

During  the  development  of  the  model  the  network  analyzer  was  calibrated  using  factory 
standards  before  fitting  the  probe. 


3.1.3  Integral  expression  of  the  admittance  (Model  3) 


Another  equivalent  formulation  for  the  admittance  of  the  line  is  also  given  by  Marcuvitz 
as  (8) 


Y  =  G+jB 


(21) 


G  = 


,  °  (  — ^[jo(ko  1fe*bsin  e)  -  Jo(ko  Ve’asin  e)f 

ln(b(a)  Ve^  sin  0 


de 


(21a) 


B=— 


V  +  b’  -  Zabcos  e))  -  2s{2K^  e’a  sit.  (|) ) 

-2Sij2koye*bsin 


de 


(21b) 


where  k,  is  the  propagation  constant  in  free  space,  a/c,  where  c  is  the  velocity  of  light  in 
vacuum  and  o)  is  the  angular  frequency;  J,  is  the  Bessel  function  of  order  zero;  and  Si  is  the 
sine  integral. 


19 


Equations  21a  and  21b  represent  an  expression  of  the  admittance  of  a  probe  as  an  integral 
over  its  aperture.  This  approximation  was  also  studied  by  Levine  and  Papas  (9),  Xu  er<z/.(10) 
and  Misra  etal.  (1 1)  and  was  derived  by  Misra  (7)  as  Equation  (18).  Therefore,  Model  2  and  3 
are  equivalent  to  each  other.  Both  models  were  obtained' from  a  variational  expression  for  the 
input  admittance  of  a  coaxial  cable  opening  into  an  infinite  ground  plane  and  semi-infinite 
sample  dielectric.  The  variational  expression  was  written  in  terms  of  an  integral  equation  for 
the  unknown  electric  field  distribution  in  the  aperture  and  was  solved  using  the  principal  (TEM) 
mode  of  the  coaxial  guide  to  represent  the  aperture  field.  Thus  the  solutions  given  by 
Equations  8  and  21a  and  21b  do  not  take  into  account  the  excitation  of  higher-order  tiKxles  in 
the  aperture. 


Implementation  of  Model  3 


Since  the  Bessel  function  of  order  zero  and  the  sine  integral  can  be  expressed  in  series 
expansions. 


Jo(z)=S 

p=K) 


(p!)2 


(22) 


Si(z)=^ 

p=0 


(.l)Pz2p+l 

(2p+l){2p+l)! 


(23) 


Equations  (21a)  and  (21b)  become 

«• 

B=y  -LKiK2<2n-|.l)i(n)f(2n+l)e*(n-H) 
n=0" 


(24) 


G=Ki 


^  ^  e*^’^i'^f^P^(-lK22)’^(b2P-a2p)(b2fl-a2‘i)S(p+q-l)| 
p=  1  2(p!q!)2 


I 


(25) 


where: 


Ki= - 1 — — 

Zoln(b/a)i(e7 


(26) 
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K2=2r/c 


(27) 


I(n)  = 


(-1)" 


(2n+l)(2n+l) 


j  ||  2(a2+b2-2abcos  0)^"i’^d0-(a2n>i+b2"+»)22"-"iS(n)j 


(28) 


(29) 


Equations  (24)  and  (2S)  can  then  be  rearranged  to  be — 


B=S  Fb(n)£*<"-"'>f2n+i) 

n=o  (30) 


«• 

p=l  q=l 


(31) 


where: 


Fb(n)=iK,K2(2"^»I(n) 

It 


(32) 


Fg(u)= 


Ki(-lK22)'^(b2p-a2p)(b2^-a2^)S(p+q+l ) 
2(p!q!)2 


(33) 


u  =p+q-2 


(34) 


It  can  be  noted  that  arrays  F^  and  are  functions  of  a  and  b  only  and  can  be  evaluated  for  a 
given  coaxial  sensor.  The  integral  in  I(n)  was  performed  by  Simpson's  three  point  rule  with 
100  subintervals.  Moreover,  a  subroutine  was  written  to  evaluate  S(n)  analytically  so  as  to 
minimize  the  numerical  errors. 


Both  the  number  of  terms  in  arrays  and  F^  were  varied  from  7  to  1 1  allowing  the 
convergence  rate  to  be  examined.  F^(l  1)  and  Fj(l  1)  were  found  to  be  less  than  2.2x10-308 
which  is  the  smallest  non-zero  real  value  that  can  be  handled  by  the  computer.  Hence  10  terms 
were  used  in  the  calculation  of  Equations  (30  and  31).  The  calculated  values  of  B  and  G  were 
substituted  into  Equation  (21)  in  order  to  get  the  values  of  admittance. 
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As  with  Mcxlel  2,  in  a  measurement  situation  a  Newton-Raphson  technique  is  used  to 
calculate  the  complex  pcnnittivity  of  the  sample  from  the  experimentally  measured  values  of  Y 
at  each  frequency  of  interest 


3.2  Experimental  Procedures 


3.2.1  Calibration 

To  measure  the  reflection  coefficient  of  a  one  port  device  an  incident  signal  I  and  its 
reflection  from  the  unknown  interface  R  should  be  measured  separately.  The  measured 
reflection  coefficient  S  ^  is  obtained  from  their  ratio. 

In  praedee  this  simple  approach  is  insufficient  because  of  imperfecdons  in  the  network 
analyzer  hardware.  In  accordance  with  network  analyzer  theory  a  measurement  is  prone  to  3 
major  sources  of  systemadc  errors: 

•  Impedance  mismatches  within  the  system  and  in  particular  improperly  matched 
source  input  impedance  which  re-reflects  signal  back  to  the  unknown  thus  adding  to  the 
original  incident  signal. 

•  The  separation  of  incident  and  reflected  signals  cannot  be  performed  with  100% 
efficiency.  This  task  is  performed  by  means  of  directional  couplers  which  have  a  finite 
directivity. 

•  Errors  may  arise  from  differences  in  the  frequency  response  of  components  of  the 
incident  and  reflection  channels. 


These  inaccuracies  generate  3  error  terms  El,  E2,  and  E3  that  are  mathematically  related  to 
the  actual  reflection  coefficient  and  the  measured  value  Tjjj. 


rn,  =  El  +E2xra/(1  -E3xra) 


(35) 


The  purpose  of  the  calibration  is  to  determine  these  errors  at  each  frequency  of  interest 
and  to  use  these  values  to  obtain  the  actual  reflection  coefficient  of  the  unknown.  These  values 
are  found  by  measuring  the  response  of  at  least  three  independent  standards  whose 


charaaerisdcs  arc  known  at  all  frequencies.  The  most  commonly  used  standards  are  a  short,  an 
open,  and  matched  loads  to  correct  for  source  match,  directivity  and  frequency  response. 


Details  of  the  calibration  procedure  are  in  the  Instruction  manual  of  the  network  analyzer  (16). 

When  open-ended  coaxial  probes  are  used  fex*  measurements  of  the  dielectric  properties  of 
materials,  small  changes  in  relatively  high  complex  reflection  coefficients  should  be  precisely 
measured,  thus  the  exact  calibration  of  the  measuring  syston  becomes  a  foreground  problem. 
Ideally  the  calibraticxi  should  be  performed  at  the  measuring  plan^nameljj  the  functional  end  of 
the  measuring  probe.  In  this  case  factory  standards  carmot  be  used. 


There  are  essentially  two  ways  of  calibrating  such  a  system* 

•  Calibration  standards  in  the  form  of  liquids  of  well-known  dielectric  pr(q)erties  may 
'  ^  used  in  conjunction  with  a  valid  probe  admittarxrc  nxxieL 

•  Alternatively,  the  calibration  may  be  performed  using  factory  standards  before 
connecting  the  probe.  A  probe  is  then  connected  at  the  calibration  port  and  the  reference  plane 
of  the  measurement  defined  by  shorting  the  end  of  the  probe  and  adjusting  the  electrical  delay 
until  the  reflection  from  the  short  circuit  at  the  functional  end  of  the  probe  exhibits  a  constant 

1 80  ®  phase  angle.  Under  these  conditions  the  network  r-  ’  ;.’zer  is  fully  calibrated  except  for 
the  reflection  from  the  conneaor  of  the  probe. 


The  advantages  of  the  first  technique  n  jst  be  weigh-xl  against  the  uncertainties  in  the 
reflection  coefficient  of  the  standards  which  depend  on  the  purity  of  the  liquid,  confidence 
interval  of  the  dielectric  parameters,  dielectric  model,  and  the  impedance  model  of  the 
functional  end  of  the  probe. 

While  the  admittance  model  is  being  assessed  it  is  not  appropriate  to  use  it  to  calibrate  the 
network  analyzer.  Therefore,  during  the  testing  of  the  riKxlels  the  latter  technique.was  adopted. 


3.2,2  Experimental  setup 


The  experimental  setup  is  shown  in  Figure  10 .  The  setup  comprises  an  HP8720  network 
analyzer  incorporating  its  own  scattering  parameters  test  set,  an  HP9000  computer  running 
Basic  5.1  and  peripheries,  a  contact  coaxial  probe  suitable  for  insertion  into  a  temperature- 
controlled  environment 
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Figure  10.  Experimental  setup  fcR*  dielectric  measurements  using  a  coaxial  probe 
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4.  RESULTS 


4.1  Comparative  Study  of  the  Three  Admittance  Models 

The  admittance  of  the  probe  was  calculated  at  a  frequency  range  ftom  130  MHz  to  20  GHz 
using  the  3  models  introduced  in  Section  3.  The  results  are  given  in  Figures  11  and  12  for  a  50 
Q  Teflon-filled  coaxial  probe  (3.5  mm  OD,  a  =  0.456  mm,  b  =  1.490  mm,  =  2.1)  in  air  and 
water  respectively.  The  complex  permittivity  of  water  was  generated  by  the  Debye  dispersion 
equation  (13,14). 


e*  =  eon  +  (ej  -0/(1  +  jtot)  (36) 

where  e_  and  e,  are  the  high  and  low  frequency  limits  of  the  relative  permittivity;  x  is  the 
relaxtition  time. 

In  Figures  11  and  12  the  squares  (□)  represent  calculations  using  Model  1,  the  stars  (*) 
Model  2  and  the  crosses  (+)  Model  3.  With  the  probe  in  air  (Fig.  11)  there  is  a  fairly  good 
agreement  between  aU  3  trodels  for  frequencies  up  to  10  GHz.  At  higher  frequencies.  Model  1 
deviates  from  the  predictions  of  the  other  models.  With  the  probe  in  water,  the  agreement 
between  the  three  models  extends  only  to  3  GHz,  Model  2  and  3  agree  up  to  7  GHz.  At  higher 
frequencies,  both  the  real  and  imaginary  parts  of  admittance  calculated  from  Model  2  alternate 
between  positive  and  negative  values  (Rg.  12).  Model  2  fails  to  converge  under  these 
conditions. 

The  admittance  of  the  probe  immersed  in  water  was  also  measured  using  the  experimental 
setup  described  in  Section  3.2.  For  the  purpose  of  these  measurements  the  network  analyzer 
was  calibrated  before  fixing  the  probe  using  factory  standards.  Figure  13  compares  the 
measured  values  with  the  predictions  of  Model  3. 
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B.  Imaginary  pan 


Figure  12 .  A  Tcflon-fiUcd  coaxial  probe  in  water. 
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A.  Measured  and  calculated  real  part. 


B.  Measured  and  calculated  imaginary  pan. 


Figure  13.  Comparison  of  the  measured  values. 


The  measurements  in  FIgtirr  13  incorporate  a  systenoatic  enor  due  to  muldple  reflecnons  by 
the  uncalibrated  probe  connector  the  effect  of  which  are  to  make  the  measured  reflection 
coefficient  to  oscillate  about  the  actual  value;  this  effea  is  mote  apparent  at  high  frequencies. 

The  agreement  between  experimental  and  calculated  admittance  values  (Fig.  13 )  is  thus  quite 
acceptable. 
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4.2  Results:  Model  1 


Using  Model  1,  good  results  are  expected  when  materials  of  low  permittivity  are  measured 
with  air  as  reference.  Equally  good  results  are  expected  when  aqueous  systems  are  measured 
using  water  as  reference.  To  test  the  performance  of  this  model  to  the  limi^a  standard  material 
of  intermediate  permitdvity  has  been  selected  for  our  study. 

Ethanediol  has  a  static  permittivity  of  ~42  at20°C  which  is  intermediate  between  the  static 
permittivity  of  water  (80. 1)  and  that  of  air.  Moreover, the  dielectric  properties  of  ethanediol  are 
well  characterised  as  a  function  of  frequency  and  temperature .  Its  dielectric  dispersion  is 
simulated  by  a  Cole-Davidson  model 


e«,  +  (£s  *0/(1 


(37) 


At  20  ®C  e5.the  static  permittivity,  is  41.4,  e„„,  the  infinite  permittivity,  is  3.7,  x,  the 
relaxation  time,  is  164.S  ps  and  P,  the  distribution  parameter,  is  0.8. 


Fig.  14  shows  the  complex  permittivity  of  ethanediol  measured  with  air  as  a  reference. 
The  solid  line  represents  the  literature  values. 
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Figure  14.  Dielectric  properties  of  ethanediol  at  20  °C. 


4.3  Results:  Model  2  and  3 


Model  2  offers  no  advantages  over  either  I  of  the  2  other  models.  Model  2  derivation  was  a 
useful  comparative  excercise  that  helped  demonstrate  the  suitability  of  Model  3.  Measurements 
reported  in  this  section  (Fig.  15  to  19)  were  obtained  using  this  model  on  ethanol  (17), 
methanol  (18),  ethanediol  (19),  formamide  (I8)and  a  1%  potassium  chloride  solution  (20). 
These  materials  exhibit  a  dielectric  behaviour  described  by  the  Debye  model  (Eq.  36),  the  Cole- 
Davidson  variation  (Equation  37)  or  by  the  Cole-(3ole  nxxlel: 

e*  =  +  (Ej  -  O  /  (1  +  (jtot)(*'“^  (38) 

An  additional  conductivity  term,  (ja/coe^p  where  is  the  permittivity  of  free  space,  is 

required  when  analysing  the  behaviour  of  materials  with  significant  ionic  conductivity.  In 
Figures  15  to  19  the  solid  line  represents  the  best  fit  of  the  data  to  the  t^)propriate  dielectric 
model  (Equations  36  to  38).  The  parameters  are  given,  with  the  95%  confidence  deviation  in 
parentheses,  they  all  fall  well  within  the  values  quoted  in  the  literamre. 
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Figure  15 .  Dielectric  properties  of  Methanol  at  22  °C.  The  solid  line  represents  Eq .  (38) 
Ej  =  34.6(0.1).  =  4.9(0.2),  x=  15.7(0.06)  ps  and  a  =  4.60.036(0.003). 


Figure  16 .  Dielectric  properties  of  Ethanol  at  22  ‘t.  The  solid  line  represents  Eq.  (36)  with 
Eg  =  25.1(0.2),  =  4.1(0.1).x=  152(0.03)  ps. 
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Figure  17.  Dielectric  properties  of  Ethanediol  at  22  ®C.  The  solid  line  represents  Eq.(37) 
Ej  =41.6(0.1,  =  3.8{0.2),  T=  134  ps  and  P  =  0.082(0.01). 


FRCOOE>JCY  ✓  H* 


Figure  18.  Dielectric  properties  of  Formamide  at  22  ®C.  The  solid  line  represents  Eq.G6) 
Ej  =109.5(0.2),  =  6.9(0.02),  t  =  37.4(0.01)  ps  and  a  =  0.052(0.003)  S/m. 


Figure  19.  Dielectric  properties  of  1%  KQ  solution  at  22  ®C.  The  solid  line  represents 
Eq.G6)  with  e,  =77.7(0.2),  =  3.7(1.43),  t  =  8.14(0.02)  ps 

and  a  '^..64(0.004). 


4.4  Error  Analysis 


The  errors  associated  with  dielectric  measuring  techniques  fall  into  2  main  categories:  (1) 
hardware  errors,  and  (2)  errors  originating  from  the  admittance  irxxlel.  The  strength  of  our 
approach  is  that  these  2  types  of  errors  were  dealt  with  independently.  For  instance,  hardware 
errors  were  corrected  using  the  manufacturers  calibration  technique  using  standards  unrelated 
to  the  admittance  model.  The  most  appropriate  admittance  model  for  the  probes  and  frequency 
range  covered  was  objectively  selected.  Then  the  measuring  technique  was  greatly  improved 
by  combining  this  model  into  the  hardware  calilsation.  The  results  (Figs.  15-19)  proclaim  the 
success  of  this  approach  by  showing  measurements  performed  on  liquids  of  known  dielectric 
properties.  In  these  results  the  measured  values  are  mostly  within  1%  of  the  calculated  values. 

We  note,  however,  that  the  uncertainties  in  e*  depend  not  only  on  the  uncertainties  in  T  at  the 

measurement  frequency,  but  also  on  the  partial  derivatives  of  E*  with  respect  to  T  which  in 

turn  depend  on  the  value  of  F  and  hence  on  the  probe  and  sample  properties.  This  complex 

situation  is  illustrated  by  plotting  F  at  a  given  frequency,  in  the  complex  permittivity  plane  (Fig. 
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20).  Regions  of  closely  packed  contours  ccHTcspond  to  large  partial  derivatives  and,  therefore, 
good  measurements  resolution.  The  use  of  such  mappings  serves  to  assess  the  measurement 
problem  and  provides  criteria  for  the  selection  of  a  probe  for  the  sample  permittivity  and 
fiequency  range  required. 


(1^) 


Figure  20.  Mapping  of  G  in  the  e*  plane  (Model  3,  3  mm  probe  at  20  GHz). 
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5.  CONCLUSIONS 


The  wcffk  reported  in  Section  2  shows  that  the  finite  element  technique  is  well  suited  to  the 
study  of  the  fringing  fields  of  coaxial  probes.  Although  a  similar  approach  has  previously  been 
used  to  calculate  the  end  of  line  c^acitance,  this  is  the  first  time  it  is  used  to  map  equipotendal 
lines.  This  approach  was  made  possible  by  the  incorporation  of  a  newly  developed  automatic 
mesh  generator  fcff  ease  of  applicadon  and  high  resoludon. 

The  use  of  a  coaxial  probe  for  dielectric  measurentents  requires  the  development  of  an  end- 
of'line  admittance  model.  A  3  mm,  polytenafluoroethylene  (PTFE)  filled  probe  was  used  in 
this  study.  The  advantages  and  limitadons  of  the  simple  lumped  parameter  model  were 
illustrated.  By  comparing  the  2  other  but  equivalent  formuladons  for  the  admittance  of  the 
coaxial  cable.  Model  2  was  found  to  fail  at  high  frequencies  (>7  GHz).  This  failure  could  be 
due  to  the  nature  of  the  exponendal  term  in  Equadon  16.  As  frequency  increases,  the  value 
of  k  increases  and  hence  the  frequency  of  vibradons  of  both  the  real  and  imaginary  parts  of  this 
term  becomes  higher.  When  the  frequency  is  increased  beyond  a  certain  limit,  the  vibradon  is 
so  rapid  that  it  becomes  meaningless  to  do  the  integradtMi  on  this  term.  Model  3  was  found  to 
be  a  better  approximadon  for  admittance  of  the  coaxial  probe  in  contact  with  lossy  media  for 
the  frequency  range  130  MHz  to  20  GHz. 


The  work  covered  by  our  report  forms  the  first  stage  in  an  in-depth  sntdy  of  dielectric 
measurements  using  contact  coaxial  probes.  This  work  inspires  further  study  involving: 

1.  Coaxial  probe  field  penetradon  in  lossy  media; 

2.  Measurements  on  finite  thickness  samples; 

3.  Measurements  on  layered  materials;  and 

4.  Developments  of  admittance  models  incorporating  the  excitadon  of  TEM  and 
higher  order  modes  at  the  probe-sample  interface. 
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